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WELL LOGGING 



2340944 



BACKGROUND OF THE ] 
Field pffrgfovgntign 

5 The invention relates to well logging, and in particular to electromagnetic 

induction well logging for determining the electrical resistivity of earth formations 

penetrated by a wellbore. More specifically, the invention is related to methods for 

* 

processing induction voltage measurements into representative values of the 
resistivity of the formations. 

10 Description of the Related Art 

Electromagnetic induction resistivity well logging instruments are well 
known in the art. Electromagnetic induction resistivity well logging instruments 
are used to determine the electrical conductivity (and its converse, resistivity) of 
earth formations penetrated by a wellbore. Measurements of the electrical 

1 5 conductivity are used for, among other things, inferring the fluid content of the 
earth formations. Typically, lower conductivity (higher resistivity) is associated 
with hydrocarbon-bearing earth formations. 

The physical principles of electromagnetic induction resistivity well 
logging are described, for example, in, H. G. Doll, Introduction to Induction 

20 Logging and Application to Logging of Wells Drilled with Oil Based Mud. Journal 

of Petroleum Technology, vol. I, p. 148, Society of Petroleum Engineers, 
Richardson TX (1949). Many improvements and modifications to electromagnetic 
induction resistivity instruments have been devised since publication of the Doll 
reference, supra. Examples of such modifications and improvements can be found, 

25 for example, in U. S. patent no. 4,837,517, U. S. patent no. 5,157,605 issued to 

Chandler et al, and U. S. patent no. 5,452,762 issued to Beard et al. 

A limitation to the electromagnetic induction resistivity well logging 
instruments known in the art is that they typically include transmitter coils and 
receiver coils wound so that the magnetic moments of these coils are substantially 

3 0 parallel only to the axis of the instrument. Eddy currents are induced in the earth 
formations from the magnetic field generated by the transmitter coil, and in the 
induction instruments known in the art these eddy currents tend to flow in ground 
loops which are substantially perpendicular to the axis of the instrument. Voltages 
are then induced in the receiver coils related to the magnitude of the eddy currents. 

3 5 Certain earth formations, however, consist of thin layers of electrically conductive 
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maierials interleaved with thin layers of substantially non-conductive material. 
The response of the typical electromagnetic induction resistivity well logging 
instrument will be largely dependent on the conductivity of the conductive layers 
when the layers are substantially parallel to the flow path of the eddy currents. The 
5 substantially non-conductive layers will contribute only a small amount to the 
overall response of the instrument and therefore their presence will typically be 
masked by the presence of the conductive layers. The non-conductive layers, 
however, are the ones which are typically hydrocarbon-bearing and are of the most 
interest to the instrument user. Some earth formations which might be of 
It) commercial interest therefore may be overlooked by interpreting a well log made 
using the electromagnetic induction resistivity well logging instruments known in 
the art. 

One solution to the limitation of the induction instruments known in the art 
is to include a transverse transmitter coil and a transverse receiver coil on the 

1 5 induction instrument, whereby the magnetic moments of these transverse coils is 
substantially perpendicular to the axis of the instrument. Such a solution was 
suggested in, L. A. Tabarovsky and M. I. Epov, Geometric an d Frequency 
Focusing in Ex ploration of Anisotropic Seams, Nauka, USSR Academy of 
Science, Siberian Division, Novosibirsk, pp. 67-129 (1972). Tabarovsky and Epov 

2 0 suggest various arrangements of transverse transmitter coils and transverse receiver 
coils, and present simulations of the responses of these transverse coil systems 
configured as shown therein. Tabarovsky and Epov also describe a method of 
substantially reducing the effect on the voltage induced in transverse receiver coils 
which would be caused by eddy currents flowing in the weilbore. The wellbore is 

2 5 typically filled with a conductive fluid known as drilling mud. Eddy currents 

which flow in the drilling mud can substantially affect the magnitude of voltages 
induced in the transverse receiver coils. The wellbore signal reduction method 
described by Tabarovsky and Epov can be described as "frequency focusing", 
whereby induction voltage measurements are made at more than one frequency, 

3 o and the signals induced in the transverse receiver coils are combined in a manner 

so that the effects of eddy currents flowing within certain geometries, such as the 
wellbore, can be substantially eliminated from the final result. Tabarovsky and 
Epov, however, do not suggest any configuration of signal processing circuitry 
which could perform the frequency focusing method suggested in their paper. 
35 A device which can measure "frequency focused" transverse induction 
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measurements is described in co-pending patent application serial no. 08/686,848 
filed on July 26, 1996, entitled, "Method and Apparatus for Transverse 
Electromagnetic Induction Logging", and assigned to the assignee of this 
invention. Interpreting measurements made from conventional induction logging 
5 instruments such as described in U. S. patent no. 4,837,517, U. S. patent no. 

5,157,605 issued to Chandler et al, and U. S. patent no. 5,452,762 issued to Beard 
et al typically involves a process known as inversion. Inversion includes 
generating an initial estimate of the probable spatial distributions of resistivity 
around the logging instrument, and using the estimated spatial distribution to 

ID generate an expected response of the particular logging instrument given the 
estimated spatial distribution of resistivity. Differences between the expected 
response and the measured response are used to adjust the model of spatial 
distribution. The adjusted model of spatial distribution is then used to generate a 
new expected instrument response. The new expected response is then compared 

15 to the measured response. This process is repeated until the difference between the 
expected response and the measured response reaches a minimum. The apparent 
spatial distribution of resistivity which generates this "closest" expected response 
is deemed to be the distribution which most probably represents the spatial 
distribution of resistivities in the earth formations surveyed by the induction 

2 0 logging instrument. See for example U. S. patent no. 5,703,773 issued to 
Tabarovsky et al. 

A limitation to using inversion methods for processing signals such as from 
the instrument described in patent application no. 08/686,848 is that the signals 
represent such complex phenomena as to make an initial estimate of the spatial 

2 5 distribution of resistivity impracticable to determine. * 

SUMMARY OF THE INVENTION 
Embodiments of the invention can provide methods for generating an 
improved estimate of horizontal conductivity, dip angle, azimuth and anisotropy 

3 0 parameter of an earth formation penetrated by a wellbore. The improved estimate is 

generated from dual-frequency transverse electromagnetic induction measurements. 
The method comprises generating an initial estimate of the horizontal conductivity, 
dip angle, azimuth and anisotropy parameter from the dual-frequency transverse 
induction measurements made at each one of a plurality of base frequencies. The 
3 5 initial estimates from each of the plurality of base frequencies are input into a primary 
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trained neural network. The improved estimate ofhorizontal conductivity, anisotropy 
parameter, dip and azimuth is calculated by the trained neural network. 

The neural network is trained by generating models of earth formations, each 
having a known value ofhorizontal conductivity, anisotropy parameter, dip angle and 
5 azimuth. Voltages which would be measured by the transverse electromagnetic 
induction instrument in response to earth formations having the physical parameters 
of each model are synthesized Initial estimates from the synthesized voltages are 
calculated, .and the initial estimates and known values from each of the models are 
input to the neural network to cause it to leam a relationship between the initial 
10 estimates and the known values in the models. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 shows an induction instrument disposed in a wellbore penetrating 
earth formations. 

1 5 Figure 2 shows a functional block diagram of the induction instrument of the 

invention. 

Figure 3 A shows the transmitter coil portion of the coil mandrel unit of the 
instrument in more detail 

Figure 3B shows the receiver coil portion of the coil mandrel unit of the 
2 0 instrument in more detail 

Figure 4 shows a functional block diagram of a transmitter controller and 
signal generator for the instrument. 

Figure 5A shows a graph of the component frequencies of the transmitter 
current. 

2 5 Figure 5B shows a graph of the composite waveform of the transmitter current. 

Figure 6A show a graph of the voltage induced in the receiver coil as a result 
of the current shown in Figure 5B flowing through the transmitter coil. 

Figure 6B shows the components of the voltage induced in the receiver and 
how digital samples made at certain times represents the difference in peak amplitude 

3 o between the two components of the induced voltage. 

Figure 7 shows the logging instrument inserted in a wellbore to explain the 
relative orientation of the instrument with respect to anisotropy axes of the earth 
formations. 

Figure 8 is a flow chart of the process of generating initial estimates of the 
3 5 conductivity distribution surrounding the logging instrument. 
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Figure 9 shows a flow chart of the process of generating improved estimates 
using the trained neural network. 

Figure 10 shows a flow chart of the process of training the neural network. 

5 DESCRIPTION OF THE PREFERRED EMBODIMENT 

L Apparatus for measuring transverse induc tion signals 
Figure 1 shows an electromagnetic induction well logging instrument 10 
disposed in a wellbore 2 drilled through earth formations. The earth formations are 
shown generally at 4. The instrument 1 0 can be lowered into and withdrawn from the 
10 wellbore 2 by means of an armored electrical cable 6 or similar conveyance known 
in the art. The instrument 10 can be assembled from three subsections: an auxiliary 
electronics unit 14 disposed at one end of the instrument 10; a coil mandrel unit 8 
attached to the auxiliary electronics unit 14; and a receiver/signal processing/telemetry 
electronics unit 12 attached to the other end of the coil mandrel unit 8, this unit 12 
1 5 typically being attached to the cable 6. 

The coil mandrel unit 8 includes induction transmitter and receiver coils, as 
will be further explained, for inducing electromagnetic fields in the earth formations 
4 and for receiving voltage signals induced by eddy currents flowing in the earth 
formations 4 as a result of the electromagnetic fields induced therein. 

2 0 The auxiliary electronics unit 14, as will be further explained, can include a 

signal generator and power amplifiers to cause alternating currents of selected 
frequencies to flow through transmitter coils in the coil mandrel unit 8. 

The receiver/signal processing/telemetry electronics unit 12, as will be further 
explained, can include receiver circuits fordetecting voltages induced in receiver coils 

25 in the coil mandrel unit 8, and circuits for processing these received voltages into 

signals representative of the conductivities of various layers, shown as 4A through 4F 
of the earth formations 4. As a matter of convenience for the system designer, the 
receiver/signal processing/telemetry electronics unit 12 can include signal telemetry 
to transmit the conductivity-related signals to the earth's surface along the cable 6 for 

30 further processing, or alteratively can store the conductivity related signals in an 
appropriate recording device (not shown) for processing after the instrument .10 is 
withdrawn from the wellbore 2. 

The electrical configuration of the instrument 10 can be better understood by 
referring to a functional block diagram of the instrument 10 shown in Figure 2. The 

3 5 auxiliary electronics unit 14 can include a transmitter controller 24 and a combination 
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analog to digital converter/digital signal processing unit (ADC/DSP) 26, both of 
which are preferably enclosed in a thermal insulating flask 28. The flask 28 can be 
of a type known in the art and is provided to maintain stable temperature, and 
consequently stable frequency, of the transmitter controller 24 and ADC/DSP 26. The 
transmitter controller 24 and ADC/DSP 26 preferably receive electrical power from 
a DC-DC converter 30. The electrical power is preferably conducted along a power 
line 30A as direct current, so that as the power passes through the mandrel unit 8 on 
the way to the auxiliary electronics unit 14 the electrical power will not materially 
increase the amount of stray voltage induced in coils in the mandrel unit 8. The 
transmitter controller 24 can include a signal generator, which will be further 
explained, for generating an alternating voltage signal at two different frequencies. 
An analog signal output 24A of the transmitter controller 24 can be connected to a 
timer controller 22 which selectively operates, at an appropriate time as will be further 
explained, each of three power amplifiers 1 6, 1 8, 20. The output of each of the power 
amplifiers 1 6, 1 8, 20 is connected to one corresponding transmitter coil set (not shown 
in Figure 2) in the mandrel unit 8. 

The ADC/DSP 26 can be connected to a reference tap on the output of each 
power amplifier 1 6, 1 8, 20. A portion of the current flowing through each transmitter 
coil (located in the mandrel-unit 8) from power amplifiers 16, 18, 20 is conducted 
provide a transmitter current reference for the transmitter controller 24, and for 
receiver circuits located in the receiver/signal processing/telemetry electronics unit 12, 
as will be further explained. The current so detected can be digitized in the ADC/DSP 
26 to provide the transmitter current reference in digital form to the transmitter 
controller 24. The use of the transmitter current reference will be further explained. 

The receiver/signal processing/telemetry electronics unit 12 can include 
preamplifiers 32, 34, 36, 38, 40 each of which is connected to one of the receiver coil 
sets (which will be further explained) in the coil mandrel unit 8. The output of each 
preamplifier can be connected to a corresponding analog-to-digital converter/digital 
signal processor (ADC/DSP), shown as 54, 52, 50, 48, 46 wherein the output of each 
preamplifier 32, 34, 36, 38, 40 is digitized and processed into a signal corresponding 
to the voltages induced in the corresponding receiver coil (not shown in Figure 2) to 
which each preamplifier is connected. Timing of operation for the ADC/DSP circuits 
54, 52, 50, 48, 46 can be provided by a controller 56. Preferably, controller 56 
operates the ADC/DSP circuits 54, 52, 50, 48, 46 so that digital signal samples are 
made by the ADC portion of each ADC/DSP circuit at a predetermined time with 
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respect to the generation of the alternating current flowing through the transmitter 
coils. The time can be determined by a clock and synchronization signals conducted 
over control line 30B from the transmitter controller 24. The controller 56 preferably 
times digitization from each ADC/DSP circuit so that the digital samples are 
5 synchronized with respect to the same signal phase in each cycle of the alternating 
voltage induced in each receiver coil. In this manner, the signal samples can be 
synchronously stacked to reduce noise in the signal output from each ADC/DSP 
circuit. A method of synchronous stacking signal digital signal samples to reduce 
noise is described in U. S. patent no. 5,452,762 issued to Beard et al. The ADC/DSP 

10 circuits 54, 52, 50, 48, 46 in the receiver/signal processing/telemetry electronics unit 
12 can be similar in design to the ADC/DSP 26 in the auxiliary electronics unit 8 as 
a matter of convenience for the system designer. 

The receiver/signal processing/telemetry electronics unit 12 can also include 
a calibration circuit 42 and an associated ADC/DSP circuit 44 connected thereto. A 

15 portion ofthe alternating current signal used to drive the power amplifiers 16, 18,20 
can be conducted to the calibration circuit 42 over analog signal line 30C. Analog 
signal line 30C is preferably electrostatically shielded to reduce parasitic induction of 
the alternating current signal into the receiver coils in the coil mandrel unit 8. On 
command from the controller 56, the calibration circuit 42 periodically conducts a 

2 0 sample ofthe alternating current to each ofthe receiver preamplifiers 32, 34, 36, 38, 
40. Since the alternating current signal thus conducted to the preamplifiers is in each 
case substantially identical, small differences in response characteristics of each 
preamplifier can be determined. The alternating current signal conducted to the 
preamplifiers is also digitized in a separate ADC/DSP 44 to generate a reference 

2 5 signal for determining the response characteristics of each preamplifier. The digitized 

output of each preamplifier from ADC/DSP's 46-54 is conducted, along with the 
digitized reference to the controller, where the response of each preamplifier can be 
determined as the change in the reference signal corresponding to each preamplifier 
when compared to the reference signal. Any necessary adjustments to the response 

3 0 ofthe preamplifiers 46-54 may be performed numerically by adjusting the acquisition 

tuning and numerical gain applied to digital samples from each ADC/DSP to match 
the measured difference in response between the reference signal and the output of 
each ofthe preamplifiers 46-54. This response calibration system is provided so that 
the measurements of the voltages induced in each receiver coil will be less affected 
35 by variations in response of each of the preamplifiers. 
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The controller 56 receives digital signal samples from each ADC/DSP 
connected to it and calculates the magnitudes of the voltages induced in each one of 
the receiver coils in the mandrel unit 8 based on the output of the respectively 
interconnected ADC/DSP's 54, 52, 50, 48, 46, 44. The induced voltage magnitudes 
5 thus calculated in the controller 56 may be conducted to a telemetry interface 58 for 
insertion into a signal telemetry format provided by a telemetry transceiver 60. The 
telemetry transceiver 60 can transmit signals to the earth's surface corresponding to 
the .calculated magnitudes. Alternatively, magnitude values calculated in the 
controller 58 may be stored in an appropriate recording device (not shown) for 

10 processing after the instrument 10 is withdrawn from the wellbore (2 in Figure 1). 

The arrangement of transmitter coils and receiver coils on the coil mandrel unit 
8 can be better understood by referring to Figures 3A and 3B. The transmitter coil 
section of the coil mandrel unit 8 is shown in Figure 3 A. A transmitter coil which can 
be wound so that its axis, and thereby its magnetic moment, is along an axis XI is 

1 5 shown at TX. Axis XI by convention will be referred to as being parallel to the X- 

axis. Coil TX is preferably substantially perpendicular to the axis of the instrument 
(10 in Figure 1). Coil TX can be electrically connected to the output of one of the 
power amplifiers (such as 16 in Figure 2). When alternating current flows through 
transmitter coil TX, an alternating electromagnetic field is in induced, which causes 

2 0 eddy currents to flow in "ground loops" in the wellbore (2 in Figure 1 ) and in the earth 

formation (4 in Figure 1) substantially coaxially about axis XI and parallel to the axis 
of the mandrel unit 8 and the instrument (10 in Figure 1). 

A short distance along the axis of the coil mandrel unit 8 can be another 
transmitter coil TZ. Coil TZ can be wound so that it axis Zl is substantially parallel 
25 to the axis of the instrument 1 0 (which by convention is generally referred to as the 
Z-axis). Coil TZ can be connected to the output of another one of the power 
amplifiers (such as 20 in Figure 2). Alternating current passing through coil TZ 
induces eddy currents in the wellbore 2 and formation 4 which flow in ground loops 
substantially coaxial with axis Zl and substantially perpendicular to the axis of the 

3 0 mandrel unit 8. 

Located a short distance further along the axis of the mandrel unit 8 can be a 
mutual balancing or "bucking" coil BX, corresponding to the X-axis transmitter coil 
TX. The winding axis X2, and therefore the magnetic moment, of coil BX can be 
substantially parallel to the axis XI of coil TX. Coil BX can be series connected in 
3 5 opposite polarity to coil TX, between coil TX and power amplifier 1 6. Bucking coil 
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TX provides that the output of a corresponding X-axis receiver coil (which will be 
further explained) is substantially zero when the instrument is disposed in a non- 
conductive medium such as air. As is understood by those skilled in the art, using 
"bucking" coils to null the corresponding receiver coil output in a non-conductive 
5 environment can be performed either by providing such bucking coils connected in 
series with the corresponding receiver coil, or alternatively can be connected in series 
with the transmitter coil In the present embodiment of the invention it is preferable 
to provide a bucking coil in series with the corresponding transmitter coil to simplify 
impedance matching between the corresponding receiver coil and its associated 

10 preamplifier (such as 44 in Figure 2), and thereby to improve the ability of the 
circuitry associated with each receiver coil to handle signals over a wide frequency 
range. The reactance of a bucking coil and its associated wiring would complicate 
impedance matching, and adjusting for signal response characteristics for a wide band 
response receiver coil, when connected in series with the receiver coil because this 

1 5 reactance is frequently nearly the same as the reactance of the receiver coil 

Still another short distance along the axis of the mandrel unit 8 is a F-axis 
transmitter coil TY. Coil TY is preferably wound so that its axis Yl, and therefore 
its magnetic moment, are substantially perpendicular to both the axis of the instrument 
10 and to the magnetic moinent of coil TX. Coil TY can be connected to power 

2 0 amplifier 1 8. Alternating current flowing through coil TY induces a magnetic field 

which causes eddy currents to flow in the wellbore 2 and the earth formation 4 in 
ground loops substantially coaxial with axis Yl and parallel to the axis of the 
instrument 10. The eddy current ground loops.corresponding to coil TY would also 
be substantially perpendicular to the ground loops associated with coils TX and TZ 
25 if the coils TX, TY, TZ are arranged as described herein. 

Bucking coils associated with transmitter coils TZ and TY are shown at BZ 
and BY, respectively. Bucking coils BZ and BY are electrically connected between 
their respective transmitter coils TZ, TY and power amplifiers 20, 18 in opposite 
polarity, as is bucking coil BX. Bucking coil BZ is wound to have its axis and 

3 0 magnetic moment along Z2 and BY is wound to have its axis and magnetic moment 

along Y2. Z2 is substantially parallel to Zl, and Y2 is substantially parallel to Yl. 

A suitable arrangement of receiver coils for the invention is shown in Figure 
3B. At the lowermost end of the receiver coil section of the coil mandrel unit 8 can 
be an .XT-axis receiver coil RX. Coil RX can be wound so that its sensitive direction 
35 is parallel to axis XI as for transmitter coil TX (shown in Figure 3 A). Eddy currents 
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flowing in ground loops corresponding to coil TX will induce voltages in coil RX 
proportional in magnitude to the magnitude of the previously explained TX-related 
eddy currents. The eddy currents themselves are proportional to the electrical 
conductivity in the path of these ground loops. 

A short distance along the axis of the coil mandrel unit 8 is a Z-axis receiver 
coil RZ wound to have its sensitive direction substantially parallel to Zl, as for its 
corresponding transmitter TZ. Eddy currents flowing in the previously explained 
ground loops related to coil TZ will induce voltages in coil RZ proportional to the 
magnitude of these eddy currents. 

The mandrel unit 8 can include a F-axis receiver coil having a sensitive 
direction parallel to Yl and is shown at RY. Eddy currents associated with coil TY 
will induce similar type voltages in coil RY. 

If the layers of the earth formations (4A through 4F in Figure 1) are 
substantially perpendicular to the axis of the mandrel unit 8, then measurements made 
by the J?-axis coils in combination with measurements made by either the X- or F-axis 
coils would be sufficient to resolve anisotropy of the conductivity of the earth 
formations. It is frequently the case, however, that the layers 4A-4F are not 
perpendicular to the axis of the mandrel unit 8 either because the wellbore (2 in Figure 
1) is inclined from vertical, <5r the layers 4A-4F are not horizontal (referred to in the 
art as "dipping beds") or a combination of these two factors. Therefore in order to 
resolve the anisotropy of the conductivity, the coil mandrel unit 8 of the invention 
preferably includes cross-axial receiver coils. One such cross-axial receiver coil is 
shown at CXY. Coil CXY receives voltages, induced as a result of eddy current 
magnetic fields which are parallel to the Yl axis (parallel to the magnetic moment of 
the Y-axis transmitter coil TY). These eddy currents may be induced as a result of 
current flowing through transmitter coil TX. As previously explained, coil TX 
includes bucking coil BX to null the output of receiver coil RX in a non-conductive 
environment. Since coil CXY is located at a different axial spacing than coil RX, 
however, nulling the output of coil CXY would require a bucking transmitter coil 
located at a different axial position than coil BX. As a matter of convenience for the 
system designer, the output of coil CXY can be nulled by including a receiver bucking 
coil connected in series and opposite polarity with coil CXY. This receiver bucking 
coil is shown at BXY. Methods of adjusting the axial position of receiver bucking 
coils such as BXY to null the output of the corresponding receiver coil CXY are well 
known in the art. The present embodiment of the invention includes cross-component 
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coil CXY instead of merely using receiver coil RY for the same reason as coil CXY 
includes associated bucking coil BXY, namely that nulling the output of receiver coil 
RY to match transmitter coil TX in a non-conductive environment would require the 
use of an additional bucking coil for cross component detection, as well as the original 
5 bucking coil for direct detection of tfie signal from its associated transmitter coil (TY 
in this case). As a matter of convenience for the system designer the present 
embodiment includes separate cross-component coils such as CXY. It is to be 
understood that RY could be used for cross-component detection when combined with 
an appropriate bucking coil, and therefore the use of separate cross-component coils 

10 should not be construed as a limitation on the invention. 

Another cross-axial receiver coil which can be included in the invention is 
shown at CXZ. Coil CXZ receives voltages induced along the Z-axis caused by eddy 
currents flowing in the earth formation as a result of current flowing through theX- 
axis transmitter coil TX (along XI). Coil CXZ can include a receiver bucking coil 

1 5 BXZ similar in function to bucking coil BXY. Adjusting the combined output of coils 
BXZ and CXZ to be zero in a non-conductive environment can be performed in a 
similar manner to that used to null the combined output of coils CXY and BXY in a 
non-conductive environment. 

The electrical connections between the receiver coils and the receiver/signal 

2 0 processing/telemetry electronics unit (12 in Figure 2) can be better understood by 
referring once again to Figure 2. Receiver coil RX can be connected to the input of 
preamplifier 32. Receivers RY and RZ are connected, respectively, to the inputs of 
preamplifiers 34 and 36. Cross-axial receiver.coil CXY and bucking coil BXY are 
series connected to the input of preamplifier 38. Cross-axial receiver coil CXZ and 

2 5 bucking coil BXZ are series connected to the input of preamplifier 40. Preamplifiers 

32-40 are each selected to provide a signal output level compatible with the dynamic 
range of analog to digital converter portion of the ADC/DSP circuit connected the 
output of each preamplifier. As previously explained, each preamplifier 32, 34, 36, 
38, 40 can be connected to an associated ADC/DSP 54, 52, 50, 48, 46. 

3 0 The ADC/DSP's 54, 52, 50, 48, 46 each generate digital samples of the output 

of the preamplifier connected to it. The acquisition of the digital samples is timed by 
the controller 56. The controller 56 can be programmed to cause each associated 
ADC/DSP to generate digital samples of the output of the corresponding preamplifier. 
The controller 56 commands each ADC/DSP to generate a plurality of samples during 
35 each cycle of the alternating current flowing through each of the transmitter coils. 
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These digital signal samples can be timed to have a predetermined phase with respect 
to the alternating voltages induced in each of the receiver coils RX, RY, RZ, CXY, 
CXZ. The significance of the timing of the digitization will be further explained. 

The hardware configuration of the instrument 1 0 having been explained, the 
timing and control of the power amplifiers (16, 18, 20 in Figure 2) and the 
ADC/DSP's (54, 52, 50, 48, 46 in Figure 2) will be explained in more detail. 
Referring now to Figure 4, the transmitter controller 24 can include a read only 
memory (ROM) 62 which contains a. digital representation of the desired waveform 
of the current to be passed through" transmitter coils (TX, TY, TZ in Figure 3 A). The 
digital representation typically consists of numbers corresponding to the magnitude 
of the desired waveform sampled at spaced apart time intervals. The output of the 
ROM 62 is timed by a clock 64, which may form part of the transmitter controller 24, 
so that the numbers exit the ROM 62 at spaced apart time intervals and are conducted 
to a digital to analog converter (DAC) 66. The DAC 66 converts the numbers 
conducted sequentially from the ROM 62 into corresponding fractional amounts of 
a reference voltage source [V^ 70 connected to the DAC 66. The output of the DAC 
66 then consists of analog voltages proportional to the numbers input from the ROM 
62. Since the output of the DAC 66 changes in step with each new number conducted 
from the ROM 62, the DAG66 would appear if graphed as a series of "stair-steps" 
The output of the DAC 66 is therefore preferably conducted to a low-pass filter 68 to 
smooth the "stair-step" like output of the DAC 66 into a continuous, smooth 
waveform. The output of the filter 68 can be conducted to the input of each power 
amplifier (16, 18, 20 in Figure 2). It is to be understood that using the digital circuit 
just described herein for generating a driver signal for the power amplifiers 1 6, 1 8, 20 
is a matter of convenience for the system designer and is meant only to serve as an 
example of circuits which could generate the desired transmitter current waveform. 
Analog signal generator circuits could as easily perform the required signal generation 
function. 

As previously explained, a reference tap on each power amplifier 16, 18, 20 
conducts a portion of the transmitter current to the ADC/DSP 26 in the auxiliary 
electronics unit 14. The ADC/DSP 26 generates digital samples of the transmitter 
current and conducts the samples to the transmitter controller 28. The transmitter 
controller 28 can calculate differences between the digitized samples of the transmitter 
current and the numbers stored in the ROM 62. These differences can include 
changes in amplitude and phase of the transmitter current with respect to the desired 
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amplitude and phase of the transmitter current. These differences can be used to 
generate adjustment factors for the numbers stored in the ROM 62 so that the desired 
amplitude and phase can be more closely generated in the transmitter current. It is to 
be understood that analog circuitry known in the art can be used to perform the 
5 adjustments to the transmitter current waveform as just described. The use of the 
digital circuitry described herein for adjusting the transmitter current waveform is a 
matter of convenience for the system designer and is not meant to limit the invention. 
The changes calculated between the numbers in the ROM 62 and the digitized 
transmitter current can also include a number of cycles of the clock 64, whereby can 

10 be determined the actual phase of the transmitter current with respect to the apparent 
phase of the transmitter current waveform as synthesized by the numbers output from 
the ROM 62. It is contemplated that the clock 64 can have a sufficiently high 
frequency whereby this phase difference can be determined to a very high degree of 
precision. The number of clock cycles of phase difference can be conducted to the 

15 controller (56 in Figure 2) in the telemetry unit (12 in Figure 2) over a serial link, 
shown in Figure 2 as 30B. The clock 64 can be used to operate both the transmitter 
controller 28 and the controller (56 in Figure 2) so that generation of digital signal 
samples of the receiver voltages can be more accurately synchronized to the 
transmitter current 

20 A method of signal processing known as "frequency focusing" enables 

determination of the conductivity of the earth formations, particularly in directions 
along the X- and K-axes, while substantially excluding the effects of eddy currents 
flowing in the weilbore (2 in Figure 1). In a novel aspect of the invention, frequency 
focusing can be performed by having the transmitter current waveform include 

25 sinusoids at two different frequencies, combined in a predetermined relationship of 
amplitude and phase between each component frequency. The frequency for the 
transmitter current can be within a range of about 10-70 KHz for the first frequency 
and about 30-2 1 0 KHz for the second frequency, as will be further explained. 

In the present embodiment of the invention, the transmitter current waveform, 

3 0 represented by J p should follow the relationship: 

h = *T1 + *T2 CO 

where l„ and / n represent, respectively, the transmitter current waveforms at the first 
o„ and the second <u, component frequencies, and where the relative amplitudes of 
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/„ and I„ follow the relationship: 



( (2) 
— I Woy) 



(0- 



where I 0 represents an arbitrary reference magnitude, typically proportional to the 
levsl of V ref (70 in Figure 4). It is desirable for to be an odd-number harmonic 
multiple of cj Jy and in the" preferred embodiment, 0% is the third harmonic. The 
5 transmitter current waveforms at the two component frequencies should have the same 
initial phase (zero) at the beginning of each cycle of the transmitter current at the first 
frequency. It has been determined that if the transmitter current follows the 
magnitude relationship described in equation (2), then the desired signal 
characteristics of the voltages induced in the receiver coils (RX, RY, RZ, CXY, CXZ 

10 in Figure 3B) can be determined by directly measuring components of the induced 
voltages which have a certain time relationship with respect to the cuiTent flowing 
through the transmitter coils (TX, TY, TZ in Figure 3 A) at the first frequency. The 
components, at these times, of the voltages induced in the receiver coils by a two- 
frequency magnetic field having the frequency, phase and amplitude relationship 

15 described in equations , and 2 are inherently substantially insensitive to voltages 
induced by eddy currents flowing in the wellbore (2 in Figure 1) and are substantially 
correspondent only to the magnitude of the eddy currents flowing only in the earth 
formations. By selecting the two component frequencies and relative amplitudes for 
the transmitter current waveform as described in equation (2), the conductivity of the 

2 0 earth formation can be directly related to the difference between the components of 
the induced voltages at each component frequency. 

In this embodiment of the invention, the difference in the magnitudes of the 
components of the induced voltages at the two frequencies can be measured directly 
by programming the controller (56 in Figure 2) to time acquisition of digital signal 

2 5 samples, represented by /„, to occur twice during each full cycle of the transmitter 
voltage at the first (lower) frequency according to the expression: 



(0, 

The digital samples of the induced voltages in the receiver coils made at these times 
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will directly represent the difference in magnitude between the components of the 
induced voltage at each frequency. The digital signal samples made at these times can 
then be directly related to the conductivity of the earth formations. 

The manner in which the magnitude of these digital signal samples directly 
5 represents the difference between the induced voltage magnitudes at the first and 
second component frequencies can be better understood by referring to Figures 5 A, 
5B, 6A and 6B. Figure 5 A shows a graph of each of the two frequency components 
of the alternating current flowing through the transmitter coil. The current magnitude 
at the first frequency is represented by curve J r/ , and the current magnitude at the 

10 second frequency is shown by curve J n . As previously explained, the second 
frequency can be the third harmonic multiple of the first frequency and have an 
amplitude relationship as previously described in equation ( 1 ). The composite current 
waveform is shown in Figure 5B as I T The voltage which is induced in the receiver 
coil as a result of eddy currents flowing in the formation is shown in Figure 6A, 

1 5 wherein the eddy currents are induced by the magnetic field generated by the current 
(J r in Figure 5B) flowing through the transmitter coil. Digital signal samples can be 
made at times shown in Figure 6A. Sample 1 is shown as timed to be at one-quarter 
cycle at the first frequency (1/2 n / <*),). This time corresponds to n = 0 in equation 
(3). Sample 2 is shown timed at three-quarter of the cycle at the first frequency (3/2 

20 n I G),), which corresponds to n = 1 in equation (3). The reason that digital samples 

made at these relative times represent the difference in magnitudes between the 
receiver voltage components at the first and at the second frequencies can be better 
understood by referring to Figure 6B, which shows the receiver voltage as its 
individual frequency components: at the, first frequency, shown by curve V R @ o,; 

25 and at the second frequency, shown by curve V R @ o> 2 . If the alternating current were 
applied at each frequency individually to the transmitter coil, the voltage induced in 
the receiver coil. would be shown by the individual component curves as in Figure 6B. 
If the current at the second frequency has the frequency and timing relationship with 
respect to the current at the first frequency, as described herein, the induced voltage 

30 at the first frequency will reach a peak value at the times at which the induced voltage 
at the second frequency will reach a peak value but at the opposite polarity. Since the 
two frequencies of current are superimposed (passed through the transmitter 
simultaneously), samples of the voltage induced in the receiver coil taken at the times 
shown in Figure 6B, such as SAMPLE 1 and SAMPLE 2, will therefore directly 

35 represent the difference between the peak magnitudes of the induced voltage 
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components at the first and at the second frequency. 

Referring once again to Figure 2, when the instrument 1 0 is first energized, the 
transmitter controller 24 begins to generate a full cycle of the transmitter voltage 
waveform. The output of the transmitter controller 24, as previously explained, is 
5 conducted to the timer controller 22. In the preferred embodiment of the invention, 
X- t Y y and Z-axis measurements can be conducted sequentially. The transmitter 
controller 24 can send a command signal to the timer controller 22 to cause it to 
actuate the particular power amplifier {such asAT-axis amplifier 1 6) whose transmitter 
coil connected thereto corresponds to the axis along which the particular measurement 

10 istobemade. It is contemplated that a sufficiently precise measurement can be made 
by operating the transmitter controller through about 1,000 transmitter voltage 
waveform cycles at the first (lower) frequency, although this number of cycles is not 
to be construed as a limitation on the invention. For example, if the instrument 10 is 
to be moved relatively slowly through the wellbore (2 in Figure 1), then a larger 

1 5 number of cycles may be useable in order to obtain higher accuracy measurements. 

After the transmitter controller has operated through about 1,000 cycles, the 
transmitter controller 28 can instruct the timer controller 22 to operate another one of 
the amplifiers, such as K-axis amplifier 18, to conduct the alternating current to its 
associated transmitter coil (TY in Figure 3 A). After about another 1,000 cycles, the 

2 0 transmitter controller 28 can instruct the timer controller to repeat the process for the 
remaining (Z-axis) power amplifier 20, and after about 1,000 cycles of alternating 
current have passed through the Z-axis transmitter coil (TZ in Figure 3 A), the entire 
process can be repeated. 

During transmission from a particular transmitter coil, the controller 56 sends 

2 5 command signals to the ADC/DSP's 46-54 which are connected to the receiver coils 
which are to be detected during operation of each particular transmitter coil. For 
example, during operation of -Y-axis transmitter coil TX, ADC/DSFs 54 (connected 
to preamplifier 32, which is connected to X-axis receiver coil RX) and 48 (connected 
to preamplifier 38, which is connected to cross-component receiver coil CXY) and 46 

30 (connected to preamplifier 40, which is connected to cross-component receiver coil 
CXZ) are instructed to digitize the output of the associated preamplifiers. The 
controller 56 can instruct the corresponding ADC/DSFs to generate digital signals 
samples at the exact times described in equation (3) with respect to the transmitter 
voltage. Alternatively, the ADC/DSP's can generate at least four samples for each 

35 cycle at the second (higher) frequency. Since the controller 56 and the transmitter 
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controller 24 can be timed by the same clock (64 in Figure 4), the exact time of 
generating the digital signal samples must be adjusted by the phase delay determined 
as previously described in the transmitter controller 24. The controller 56 can delay 
sending instruction to digitize the preamplifier output by the number of clock cycles 
5 of phase delay conducted from the transmitter controller 24 with respect to a start of 
transmitter voltage cycle command. The start of transmit cycle command can also be 
sent along the serial link 30B to indicate to the controller 56 that the transmitter 
controller 24 is initiating a transmitter voltage cycle. The digital signal samples 
acquired during the transmitter voltage cycle can be synchronously stacked, as 

10 previously described herein, and can be stored in a buffer (not shown separately) in 
the controller 56 for transmission to the surface by the telemetry transceiver 60, or can 
be retained for later processing. As previously described, the samples from each 
ADC/DSP 54-46 can be processed by a discrete Fourier transform to determine the 
magnitude of the voltage components at each frequency, or the samples made at the 

15 precise times described in equation (3) can be used to determine the difference in the 
in-phase voltage components directly. 

After the previously described number of transmitter voltage cycles (which as 
previously explained can be about 1,000), the controller 56 can send digitization 
commands to the ADC/DSP associated with the receiver coil which will receive 

20 voltages induced by its axially associated transmitter coil (such as ADC/DSP 52 
associated with coils RY and TY). The receiving and digitization process can then be 
repeated for the remaining transmitter coil and axially associated receiver coil. 

2. Processing transverse induction signals into initial estimates of the 
spatial distribution of electrical conductivity around the logging instrument 

25 The principle of the invention can be illustrated by an idealized situation 

wherein the earth formation is electrically anisotropic and the well logging instrument 
such as described herein is oriented so that the magnetic dipoles of some of its 
transmitters and receivers are parallel to the principal directions (principal axes) of the 
electrical anisotropy. If the well logging instrument described herein were so 

3 0 disposed in such an anisotropic earth, the voltage measurements which would be made 
by the various receiver coils in the instrument could be directly related to the 
conductivity values of the earth formation along each of the principal directions of the 
anisotropy. Note that in this idealized situation, there would be substantially no 
relationship between the magnitudes of voltages induced in the cross-component 

3 5 receiver coils (CXY, CXZ in Figure 3B) and the conductivity of the formation. 
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The well logging instrument (10 in Figure 1), however, will typically not be 
oriented so that the magnetic dipoles of the transmitters and receivers are parallel to 
the principal directions of anisotropy of the formations. The receiver voltage 
measurements therefore cannot be directly related to the electrical conductivity along 
the principal directions of the anisotropy. Referring now to Figure 7, as is well known 
in the art the wellbore 2 may be inclined from vertical at any particular axial position 
along its length, and/or the earth formations, such as illustrated by formation layer 4A, 
may be inclined with respect to horizontal (the inclination angle known in the art as 
"dip"). A relative angle 8 may be subtended between the axis of the instrument 10 
shown as z and a "vertical" principal axis of the layer 4A shown by z ' The vertical 
principal axis is defined for purposes of this explanation as being substantially 
perpendicular to the layering attitude or bedding planes" of the earth formation, such 
as shown for layer 4A. The electrical conductivity of the formation layer 4A 
measured in a direction parallel to the vertical principal axis can be defined as the 
"vertical" conductivity, a r The angle subtended, 0, will depend on, among other 
things, both the inclination of the formation layer 4A from horizontal (with respect to 
earth's gravity) and the inclination of the wellbore 2 from vertical (with respect to 
earth's gravity). 

Similarly as for the vertical principal axis, a "horizontal" principal axis, x 'can 
be defined as being substantially parallel to the bedding plane of the formation layer 
4A. The electrical conductivity of the formation measured parallel to the horizontal 
principal axis can be referred to as the "horizontal" conductivity and can have a value 
represented by a k . The X- or F- axes of the logging instrument 10 may subtend an 
angle of rotation, represented by <|>, with respect to the plane in which the horizontal 
principal axis of the formation layer 4A lies. For purposes of explaining this 
invention, it will be assumed that the vertical principal axis z Is vertical with respect 
to the earth's gravity and that the horizontal principal axis* Is perpendicular to earth's 
gravity. It should be understood, however, that the method of the invention is 
applicable to any other orientation of the vertical principal axis z ' and horizontal 
principal axis x 'with respect to earth's gravity. 

Angle <p can be estimated directly from some of the receiver voltage 
measurements. There are five receiver voltage measurements made by the instrument 
described herein and these can be defined in terms of the transmitter magnetic field 
orientations and the receiver dipole directions they relate to such as to represent 
the RX receiver coil measurement made from the magnetic field generated by the TX 
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transmitter coil. The other receiver measurements include if^ if^ if^ fl^ The 
five receiver voltage magnetic field representations can be represented for simplicity 
by a tensor If. Angle <|> can be calculated from the AY cross-component (if^) and 
theXXand ITmagnetic field components {0^ if^ respectively) of tensor fif by the 
expression: 

tan(2<|>) = 7 H V ■ 

- 

Once the angle <J) has been estimated, an intermediate magnetic field component 
tensor If can be determined by rotation of /^through an angle of magnitude (-<f>). 
From the intermediate magnetic field component tensor, It, the angle 0can then be 
estimated by the expression: 



10 After the angle 0 is estimated, the tensor If can be converted into the 

coordinate system of the principal axes x 'and z \ by rotating through an angle of 
magnitude (- &). The resulting tensor, referred to as H represents the magnetic field 
components that would be measured if the instrument's transmitter and receiver 
dipoles were parallel to the principal axes of the earth formation layers. Particularly, 

15 the X-axis is defined as being coplanar with <£, and the Z-axis is defined as being 
parallel to a r ; . 

The voltages induced in an X-axis receiver coil (such as RX in Figure 3B) by 
the magnetic field generated by an X-axis transmitter (such as TX in Figure 3 A) can 
be expressed in terms of this magnetic field, defined as by the following 

20 ' expression: 

■ ^ [[(sinecosd)) 2 ^ -£J£. + 7 0 -(sinesin*) 2 /,! (6) 
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Similarly, the voltages induced in a F-axis receiver by the magnetic field generated 
by a F- axis transmitter, represented by can be expressed as: 

ff = -^[[(smesmtt)) 2 ^ -*J£ B *J 0 -(smBcos*)^] (7) 
4nr 3 

Z-axis receiver voltages induced by the magnetic field from a Z-axis transmitter, H & 
can be expressed as: 

^--^Icos^&K-BJi. - (8) 
4nr* 

Cross-component voltages, as measured by a receiver coil such as CXY in Figure 
3B, and resulting from the magnetic field generated by a transmitter coil such as TX 
in Figure 3 A can be represented by the expression: 

H = \AJE U + JJ raft cos* sin* (9) 
4nr 3 

Cross-component voltages, measured by a receiver coil such as CXZ in Figure 
3B from the magnetic field generated by a transmitter coil such as TX in Figure 3 A 
can be described by the expression: 



3l 

47ur 

Certain symbols used in the foregoing expressions can be defined as follows: 



H = J5L^QsmQ*m<bA 9 E 9 (10) 
47ir 3 
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A 9 = 3 ♦ 3n ♦ u 2 



B 9 = 1 + u ♦ ii 2 



^ = (11) 



In the definitions recited in equation (1 1), the following variables can be defined: 
u=-ik h r y v=-ik v s,k k 2 =-ia>fi ^K 2 ~-i6>pa„r = (/t + ey' 2 y s = (f? + fe) 1/2 
and A-k h /k w . Also in these expressions, z " r cos ff, x = r sin Ocos ^andj = r sin 
Osin (p. The various subscripted M values represent the magnetic moments of each 
5 of the transmitters (TX, TY, TZ, in Figures 3 A, 3B) in the well logging instrument 1 0. 

The "horizontal" conductivity, a ki defined previously as the conductivity in a 
direction parallel to the bedding planes, can then be determined from a low frequency 
approximation of the signal component H a (which represents the H f u component after 
it has been "rotated" into the coordinate system of the principal axes). 
1 0 For an instrument such as the one described herein where the magnetic field 

components are measured using a so-called "three coil" system (including transmitter, 
bucking and main receiver coils) and the measurements are made at two frequencies 
as described earlier herein, the horizontal conductivity a k can be determined by the 
expression: 



& apj 



(12) 



15 where represents the imaginary part of the H u component signal, and a 

represents a ratio of spacings between the particular transmitter and each of the two 
receiver coils (the main receiver coil and the bucking coil) used to measure the 
particular component signal. The vertical conductivity <r w and the anisotropy 
parameter k can be then determined from a low frequency approximation of 

20 component signal First, the anisotropy parameter A can be calculated by the 
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expression: 



J_ m 1 
k 2 = 3 



and the vertical conductivity a y can be calculated from the expression: 



(13) 



°' = § < 14 > 

It should be noted that equation (12) can be used to calculate horizontal conductivity 
from the "rotated" H u component where a well logging instrument such as described 
5 herein which makes measurements of receiver voltage at two frequencies is used. To 
practice this invention, however, it is not necessary to use an instrument such as the 
one described herein having "bucking" coils, nor is it necessary to make induction 
voltage measurements at two different frequencies. Appendices A and B include 
expressions derived for low frequency approximations using so-called "two-coil" 

1 0 systems, and for "three-coil" systems using only a single frequency alternating current 

for making measurements. So-called "two-coil" systems include only a single 
transmitter coil and a single receiver coil oriented along one of the X 9 F, and Z-axes. 
A typical three-coil system is disclosed herein, where each of the X, F, and Z-axes 
includes a transmitter coil oriented therealong, and main receiver and "bucking 11 coil 

15 oriented along each of the X, T 9 and Z-axes. As previously explained, cross 
component receivers should be included in the logging instrument. 

Figure 8 shows a flow chart of the process of the generating the initial 
estimate. Block 1 00 shows the ifieasurements made of voltages induced in each of the 
receiver coils. The voltages, indicated as V^, V^, V^, V^, by their subscripts 

2 0 represent the voltages induced in receiver coils oriented along the axis indicated by 
the second subscript letter, these voltages being induced by the transmitter coil 
oriented along the axis indicated by the first subscript letter. To practice this 
invention, it is only necessary to measure these five voltages having the indicated 
magnetic field orientations. Block 102 shows the step of estimating the angle $ from 

25 the XY, XX and YY receiver signals. After angle <J> is estimated, the measurements 
are "rotated" through angle (-<[>) to obtain the intermediate component tensor fl*, as 
shown in block 1 04. In block 1 06, angle 6 can be estimated from the XZ, ZZ and XX 
components of the intermediate vector Ft. The estimate of angle 0 is used, as shown 
in block 108, to calculate magnetic field components in the coordinate system of the 
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anisotropy a k1 a r The horizontal conductivity is calculated as shown in block 1 10 
from the magnetic field components in the anisotropy coordinate system, using a low 
frequency approximation selected from those such as shown in Appendix B, for the 
particular instrument coil arrangement and frequency measurement system actually 
used to make the induction voltage measurements. Finally, the anisotropy coefficient 
and vertical conductivity can be calculated as shown in block 1 12. 

L Generating improved values of parameters from the initial estimate. 

The process of generating an initial estimate of horizontal conductivity, 
anisotropy parameter, apparent "dip" angle and azimuth, as previously described, can 
be repeated for each of the seven "base" frequencies of the well logging instrument 
in Figure 1. The result will be seven sets of initial estimates of horizontal 
conductivity, anisotropy parameter, dip angle and azimuth. These seven sets of initial 
estimates can then be entered into a feedforward neural network to generate a final 
("improved") estimate of the horizontal conductivity, dip angle and azimuth, 

The neural network model preferred for this invention is known as radial basis 
function (RBF) because of the relatively high training speed and accuracy. An RBF 
network includes 3 "layers" of units or neurons, the input layer, the hidden (or RBF) 
layer and the output layer. When an input vector, x, is presented to the network, the 
network responds with an output vector, j>. The network output is simply a linear sum 
of certain radially symmetric functions, x zxxdy are usually related as: 

**) -GUx-XjVop; 



where G represents a radially symmetric function, usually a Gaussian function, Xj 
represents the centroid of the /-th RBF and is the mean value of the Gaussian, a } 
represents the Gaussian width and R is the vector of J?/s. w represents a coefficient 
vector of weights for the Gaussian functions. 

Using the neural network to generate the improved estimates is shown in 
Figure 9. The voltages corresponding to the five magnetic field components for each 
of the seven base frequencies are measured, shown at box 120. The measured 
voltages are used to generate initial estimates for each base frequency, as shown at 
box 122. The initial estimates are input to the network, shown at box 124, and the 
improved estimate is the output of the network, shown at box 126. 

Training the neural network includes generating a set of model earth 
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fonnations, preferably randomly, each having known values of horizontal 
conductivity, anisotropy parameter, dip and azimuth angles. This is shown in Figure 
10 at box 130. The random set of earth formation models can then be used to 
synthesize the response of the well logging instrument. The synthetic response of the 
well logging instrument should include voltages for each of the five magnetic field 
components V„, V^, V w V„, at each of the seven base frequencies of the well 
logging instrument. This is shown at box 132 in Figure 10. The synthesis will 
comprise a set of 35 voltages (corresponding to a magnetic field magnitude) for each 
"depth" level within each one ot the random set of modeled earth formations. 
Equations which can be used to synthesize the voltage measurements for each model 
are shown in Appendices A and 8 which follow the description of this invention. 

Each of the sets of synthetic voltage measurements can be used to calculate 
initial estimates, at each of the seven base frequencies, of horizontal conductivity, 
anisotropy parameter, dip and azimuth angle in the manner previously explained 
herein. This is shown inFigure lOatbox 134. After calculating the initial estimates, 
the random set of modeled earth formations, having the known values of horizontal 
resistivity, anisotropy parameter, dip and azimuth angle can be entered into the 
feedforward neural network with the calculated initial estimates in order to train the 
neural network. As a large number of modeled earth formations and initial estimates 
are entered into the neural network, the network will "learn" the relationship between 
the "true" values ofhorizontal resistivity, anisotropy parameter, dip and azimuth angle 
from the model earth formations, and the initial estimates calculated from the 
synthetic voltage measurements. This is shown in Figure 10 at box 136. 

After such training, the : neural network can be used to make more accurate 
estimates of the horizontal resistivity, anisotropy parameter, dip angle and azimuth 
angle for earth formations in an actual wellbore for which real signal voltages are 
measured during well logging operations, where the instrument passes through 
formations penetrated by the wellbore. 

It should be understood that the invention is not limited to the radial basis 
function neural network model. The invention is also not limited to using voltage 
measurements from seven base frequencies. In practice, the method of this invention 
may use more or fewer base frequency measurements, with the understanding that 
more base frequencies may improve the accuracy of the estimates, while increasing 
the amount of training time and calculation time for the more accmate estimates 
output by the neural network. 
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It should also be understood that the improved estimate may include any one 
or all of the parameters horizontal conductivity, anisotropy parameter, apparent dip 
and azimuth. The neural network may be trained or otherwise operated to put only 
one or any ones of the parameters in the output for each calculation. It is not 
necessary, for purposes if this invention, to calculate a value of the final estimate for 
all four parameters. 

Those skilled in the art will devise other embodiments of the invention which 
do not depart from the spirit of the invention as disclosed herein. 
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CLAIMS: 

1. A method for generating an unproved estimate of at least one of 
horizontal conductivity, dip angle, azimuth and anisotropy parameter of an earth 
formation penetrated by a wellbore, from dual-frequency transverse electromagnetic 

5 induction measurements of said formations, comprising: 

generating an initial estimate of said at least one of horizontal 
conductivity, dip angle, azimuth and anisotropy parameter from said dual-frequency 
transverse induction measurements, made at each one of a plurality of base 
frequencies; 

10 imparting said initial estimates from each of said plurality of base 

frequencies into a trained neural network; and 

calculating said improved estimate using said trained neural network. 

2. The method as defined in claim 1 wherein said neural network 
1 5 comprises a radial basis function model. 

3. The method as defined in claim 1 wherein training of said neural 
network comprises: 

generating models of earth formations each said model having a known 
2 0 value of said at least one of horizontal conductivity, anisotropy parameter, dip angle 
and azimuth; 

synthesizing voltages which would be measured by said transverse 
electromagnetic induction logging instrument in response to each of said models; 

calculating initial estimates from said synthesized voltages; and 

2 5 inputting said initial estimates and said known values from each of said 

models to said neural network to cause said network to learn a relationship between 
said initial estimates from said synthesized voltages and said known values. 

4. The method as defined in claim 1 wherein generating said initial 

3 0 estimate comprises: 

measuring electromagnetic induction signals induced by induction 
transmitters oriented along three mutually orthogonal axes, one of said axes 
substantially parallel to an instrument axis, said signals measured using first receivers 
each having a magnetic moment parallel to one of said orthogonal axes and second 
3 5 receivers each having a magnetic moment perpendicular to a one of said orthogonal 
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axes which is perpendicular to said instrument axis, the method comprising; 

calculating a relative angle of rotation of said one of said orthogonal 
axes which is perpendicular to said instrument axis, with respect to a direction of said 
horizontal conductivity, from said received signals measured perpendicular to said 
instrument axis; 

calculating an intermediate measurement tensor by rotating magnitudes 
of said receiver signals through a negative of said angle of rotation; 

calculating a relative angle of inclination of one of said orthogonal axes 
parallel to said axis of said" instrument, with respect to said direction of said vertical 
conductivity, from said rotated magnitudes; 

rotating said rotated magnitudes through a negative of said angle of 

inclination; 

calculating said horizontal conductivity from said magnitudes of said 
receiver signals after said second step of rotation; and 

calculating an anisotropy parameter from said magnitudes after said 
second step of rotation and calculating said vertical conductivity from said horizontal 
conductivity and said anisotropy parameter. 

5. A method for generating an improved estimate of at least one of 
horizontal conductivity, dip angle, azimuth and anisotropy parameter of an earth 
formation penetrated by a wellbore, comprising: 

inserting a transverse electromagnetic induction logging instrument 
comprising mutually orthogonal transmitters and receivers and a dual-frequency 
alternating current source selectively correctable to each of said transmitters to induce 
alternating magnetic fields in said formations along ah axis parallel thereto, said dual 
frequency source comprising a plurality of base frequencies; 

detecting voltages at each said receiver corresponding to 
electromagnetic induction phenomena in said earth formation, said voltages including 
one component corresponding to reception along the same axis as each of said induced 
magnetic fields, and reception orthogonal to said induced magnetic fields; 

generating an initial estimate of said at least one of horizontal 
conductivity, apparent dip angle, apparent azimuth and anisotropy parameter from 
said detected voltages at each one of said plurality of base frequencies; 

imparting said initial estimates from each of said plurality of base 
frequencies into a trained neural network; and 
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calculating said improved estimate using said trained neural network. 

6. The method as defined in claim 5 wherein said neural network 
comprises a radial basis function model. 

5 7. The method as defined in claim 5 wherein training of said neural 

network comprises: 

generating models of earth formations each said model having a known 
value of said at least one of horizontal conductivity, anisotropy parameter, dip angle 
and azimuth; 

10 synthesizing voltages which would be measured by said transverse 

electromagnetic induction logging instrument in response to each of said models; 

calculating initial estimates from said synthesized voltages; and 
inputting said initial estimates and said known values from each of said 
models to said neural network to cause said network to learn a relationship between 
1 5 said initial estimates from said synthesized voltages and said known values. 

8. The method as defined in claim 5 wherein generating said initial 
estimate comprises: 

measuring electromagnetic induction signals induced by induction 
20 transmitters oriented along three mutually orthogonal axes, one of said axes 
substantially parallel to an instrument axis, said signals measured using first receivers 
each having a magnetic moment parallel to one of said orthogonal axes and second 
receivers each having a magnetic moment perpendicular to a one of said orthogonal 
axes which is perpendicular to said instrument axis, the method comprising; 
2 5 calculating a relative angle of rotation of said one of said orthogonal 

axes which is perpendicular to said instrument axis, with respect to a direction of said 
* horizontal conductivity, from said received signals measured perpendicular to said 
instrument axis; 

calculating an intermediate measurement tensor by rotating magnitudes 
30 of said receiver signals through a negative of said angle of rotation; 

calculating a relative angle of inclination of one of said orthogonal axes 
parallel to said axis of said instrument, with respect to said direction of said vertical 
conductivity, from said rotated magnitudes; 

rotating said rotated magnitudes through a negative of said angle of 
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inclination; 

calculating said horizontal conductivity from said magnitudes of said 
receiver signals after said second step of rotation; and 

calculating an anisotropy parameter from said magnitudes after said 
second step of rotation and calculating said vertical conductivity from said horizontal 
conductivity and said anisotropy parameter. 

9. A method for generating estimates of at least one parameter of an earth 
formative penetrated by a wellbore, substantially as hereinbefore described with 
reference to the accompanying drawings. 
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